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摘要 :植被 光合 呼吸 模型 (VPRM ) 关键 参数 的 确定 和 优化 是 准确 计算 生态 系统 净 CO0, 交 换 (NEE) 的 基础 。 利 用 中 国 通 量 观测 
研究 联盟 ( ChinaFLUX) 长 白山 站 温带 阔 叶 红 松林 2005 年 的 通 量 观测 资料 ,对 VPRM 的 和 个 参数 (最 大 光 能 利用 率 e, .光照 为 半 
饱和 条 件 下 光合 有 效 辐射 值 PAR。 和 呼吸 参数 (ww,B) ) 进行 优化 ,并 使 用 2006 年 的 观测 资料 对 参数 优化 前 后 的 模拟 结果 进行 评 
估 。 结 果 表 明 :参数 优化 后 ,VPRM 能 够 较 好 地 模拟 长 白山 地 区 2006 FIE E T ZR NEE BTE b. Xf 30minNEE 模拟 的 平均 误 
差 为 -1.81 umol m^ s^! ,相关 系数 为 0.72 ,模拟 NEE 平均 日 变化 的 峰值 约 为 观测 值 的 97% ,相关 系数 为 0.97。 但 在 植物 非 生长 
季 模 型 对 森林 NEE 的 模拟 效果 较 差 。 模 型 模拟 30minNEE 的 平均 误差 为 0.39 umol m? s^! ,相关 系数 仅 为 0.10 ,并 且 模 拟 低 估 
NEE 平均 日 变化 白天 吸收 峰值 约 8296 ,日 变化 模拟 值 与 观测 值 的 相关 系数 为 0.50。 通 过 分 析 不 同 天 气 个 例 ,发 现 模 型 可 以 较 
好 地 模拟 晴天 条 件 下 NEE 的 变化 ,而 对 阴雨 天 NEE 的 模拟 误差 较 大 。 该 研究 有 利于 提高 VPRM 模型 对 温带 落叶 阔 叶 林 NEE 
的 模拟 能 力 ,对 进一步 改进 区 域 陆 地 NEE 的 模拟 具有 重要 意义 。 

关键 词 .植被 光合 呼吸 模型 ;生态 系统 净 CO, 交换 ;参数 优化 ;长 白山 温带 阔 叶 红 松林 
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Abstract; Optimizing key parameters of the Vegetation Photosynthesis and Respiration Model ( VPRM ) is crucial for 
accurately calculating net ecosystem CO, exchange ( NEE). The objectives of the present study were to use data measured in 
temperate broad-leaved Korean pine forest at Changbaishan ( CBS) in 2005 to optimize four VPRM parameters, maximal 
light use -éfficiency 人 es, ) half-saturation value of photosynthetically active radiation ( PAR, ), and two respiration 
parameters (a, B) band to evaluate the simulation results using data measured in 2006. Using the optimized parameters, 
the VPRM.was able to accurately simulate variation in net ecosystem CO, exchange ( NEE) during the growing season of 
2006. For 30 min NEE simulation, the mean bias was —1.81 umol m” s™', and the correlation coefficient was 0.72. In 
regards to daily variation, the peak NEE value was underestimated by 996, and the correlation coefficient was 0.97. 
However, the VPRM was unable to accurately simulate NEE during the non-growing season. During this period, the mean 


bias for 30 min NEE simulation was 0.39 umol m” s™' , and the correlation coefficient was 0.10. In regards to daily variation 
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during the non-growing season, the peak NEE value was underestimated by 8296, and the correlation coefficient was 0.50. 
Furthermore, we also found that the VPRM is better suited for simulating NEE on sunny days than on cloudy or rainy days. 
The present study facilitates application of the VPRM to studies of the NEE of temperate broad-leaved Korean pine forest 


and has important significance for improving simulation of regional terrestrial ecosystem NEE. 


Key Words: vegetation photosynthesis and respiration model ( VPRM) ; net ecosystem CO, exchange ( NEE) ; parameter 


optimization; temperate broad-leaved Korean pine forest 


系统 碳 循环 研究 的 主要 目的 之 一 是 确定 大 气 与 陆地 生态 系统 之 间 净 CO, 交换 ( NEE) 的 量 值 , 并 得 到 处 同时 空 
尺度 下 NEE 变化 的 主要 驱动 力 |。 涡 度 相关 技术 能 够 准确 地 测量 陆地 生态 系统 与 大 气 之 间 NEE 的 量 值 ,并 
在 过 去 的 20 多 年 间 , 被 广泛 用 于 全 球 不 同 生 态 系 统 类 型 NEE 的 确定 。 但 是 ,受到 观测 高 度 、 风 向 和 下 执 
面 状况 的 影响 , 涡 度 相关 技术 实际 观测 的 区 域 范围 较 小 ,一般 在 几 公里 范围 之 肉 ,得 到 的 NEE 数据 在 较 大 的 
空间 范围 上 并 不 具 代表 性 '"。 但 是 ,卫星 遥感 观测 不 但 能 够 以 固定 的 频率 对 生态 系统 要 素 进行 采样 ,其 观测 
能 够 覆盖 区 域 乃 至 全 球 范 围 。 所 以 基于 卫星 遥感 的 植被 生产 力 模型 在 估计 区 域 不 同 后 态 系统 生产 力 方面 发 
挥 着 重要 的 作用 2 。 

近年 来 ,基于 卫星 遥感 数据 的 植被 生产 力 模型 已 经 被 用 来 模拟 区 域 的 总 初级 生产 力 (GPP) 。 其 中 ,统计 
么 感 模型 发 展 较 早 ,其 原理 是 建立 地 面 观测 的 植被 生产 力 与 还 感 获取 的 植被 指数 之 间 的 统计 关系 ,该 模型 结 
构 简 单 且 统计 关系 容易 建立 。 如 Paruelo 等 利用 地 上 植被 生 诺 力 与 归 一 化 植被 植被 指数 (NDVI) 之 间 的 函数 
关系 实现 了 对 美国 中 部 草原 GPP 的 估算 "” , Beer 等 建立 的 多 变量 高 级 统计 方法 的 植被 生产 力 估算 诊断 模 
型 实现 了 对 全 球 不 同 气 候 ( 纬 度 ) 地 区 GPP 的 估算 ' 等。 然而 该 类 模型 未 考虑 植物 的 生理 过 程 ,区 域 适 用 性 
和 经 验 性 很 强 , 在 将 其 应 用 于 研究 地 以 外 的 地 区 时 需 重新 确定 其 经 验 参数 值 。 而 基于 卫星 遥感 的 光 能 利用 
率 模型 是 利用 归 一 化 植被 指数 (NDVI) 求 取 植 被 吸收 用 于 光合 作用 的 光合 有 效 辐 射 的 比例 (FAPAR ) ,结合 3》 
能 利用 率 , 在 考虑 了 植物 的 生理 过 程 林 础 上 实现 对 CPP 的 估算 。 如 第 一 代 的 CASA 模型 ( Carnegie-Ames- 
Stanford Approach ) “中 以 及 融合 嫩 涡 度 相关 观测 数据 的 CFix 模型 (Carbon Fix) 5 , EC-LUE 模型 (Eddy 
Covariance-Light Use Efficiency) ^ , 然而 ,NDVI 对 地 面 植被 指数 饱和 、 冠 层 背 景 等 要 素 十 分 敏感 '", 这 将 
导致 上 述 模型 对 GPP 的 模拟 存在 较 大 的 不 确定 性 。 因 此 , Xiao 等 人 建立 了 VPM( Vegetation Photosynthesis 
Model) 使 用 增强 型 植被 指数 (上 EVI) 代替 NDVI, 并 结合 地 面 涡 度 相关 资料 ,在 明确 生态 系统 最 大 光 能 利用 率 的 
基础 上 ,引入 温度 ,水 分 、 及 物候 变化 对 植被 光 能 利用 率 的 影响 ,从 而 实现 对 不 同类 型 植被 生产 力 的 准确 模拟 
0021 ”例如 Xiao 等 人 使 用 VPM 实现 了 对 美国 Harvard. Forest 和 Howland Forest 站 点 森林 GPP 的 估算 和 下 
Jin 等 使 用 NPM 实现 了 对 美国 7 个 作物 站 点 GPP 的 估算 ” 。VPM 虽然 能 够 较 好 地 模拟 CPP ,但 还 不 能 直接 
对 NEE 进行 模 扎 。 因 此 ,Mahadevan 等 人 建立 的 VPRM (Vegetation Photosynthesis and Respiration Model) Æ 
VPM 的 基础 此 改进 了 GPP 与 光合 有 效 辐射 (PAR ) 之 间 的 响应 关系 ,增添 了 PAR 的 非 线 性 项 ,反应 植被 对 光 
的 响应 特征 ,同时 还 增加 了 呼吸 项 ,反应 生态 系统 呼吸 对 温度 的 响应 特征 ,这 些 改进 使 模型 在 考虑 植被 生理 特 
性 的 基础 上 实现 对 NEE 的 直接 模拟 。 

但 是 ,要 实现 VPRM 模型 的 准确 模拟 ,还 要 对 模型 中 反应 植被 生理 特性 的 参数 进行 优化 。Mahadevan 等 
利用 北美 地 区 22 个 观测 站 的 通 量 观 测 资料 对 VPRM 的 参数 进行 了 优化 ,获得 了 一 组 适合 北美 地 区 不 同 植被 
类 型 的 模型 参数 ,参数 优化 后 的 VPRM 模型 较 好 地 模拟 了 各 站 点 NEE 的 时 空 变化 '”。 而 Hilton 等 根据 65 
个 北美 通 量 观测 站 点 的 观测 ,进一步 完善 了 VPRM 参数 在 北美 地 区 的 优化 '” 。 但 是 ,由 于 不 同 区 域 气候 差异 
较 大 ,和 森 林 植 被 生理 生态 学 特征 各 异 , 因 此 这 些 参 数 在 受 季 风 影 响 的 我 国 可 能 并 不 适用 。 因 此 ,为 实现 VPRM 
模型 在 我 国 不 同 区 域 的 应 用 ,就 需要 针对 我 国 不 同 气候 区 的 植被 类 型 进行 VPRM 模型 参数 的 优化 。 刘 诚 等 
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人 做 了 相关 的 研究 ,利用 千 烟 洲 亚热带 人 工 针 叶 林 的 观测 资料 对 VPRM 参数 进行 了 拟 合 优化 荆 ] 。 但 是 ,针对 
我 国 对 气候 变化 敏感 的 温带 地 区 的 森林 生态 系统 进行 VPRM 优化 及 验证 的 研究 较 少 。 

温带 针 阔 叶 林 为 我 国 4 个 主要 森林 类 型 之 一 ,在 我 国 东北 温带 气候 区 有 着 广泛 的 分 布 。 本 研究 旨 在 利用 
中 国 通 量 观测 研究 网 络 ( ChinaFLUX ) 长 白山 通 量 观 测 站 (CBS) 的 温带 阔 叶 红 松林 通 量 观测 资料 对 VPRM 模 
型 关键 参数 进行 拟 合 优化 ,并 评价 该 模型 在 温带 针 闻 叶 混 交 林 的 模拟 效果 。 从 而 获取 一 组 适合 温带 针 间 叶 混 
交 林 的 参数 ,提高 VPRM 模型 对 温带 针 阔 叶 林 NEE 的 模拟 能 力 , 并 为 未 来 在 我 国 其 他 地 区 做 VPRM 模型 参 
数 的 优化 及 验证 提供 参考 。 


1 研究 地 区 和 研究 方法 


1.1 研究 区 域 概况 
长 白山 温带 红 松 阔 叶 林 通 量 观测 站 位 于 吉林 省 延边 朝鲜 族 自治 州 安 图 县 二 道 白河 镇 。 该 站 地 处 长 白山 
自然 保护 区 内 ,其 长 白山 站 点 地 理 .气候 和 植被 状况 如 表 1 所 示 。 


表 1 长 白山 站 点 地 理 ,气候 和 植被 状况 
Table 1 The geography, climate and vegetation condition in ChangBai station 


项 目 Projects 描述 Description 项 目 Projects 描述 Description 
地 点 Location 吉林 省 延边 朝鲜 族 自治 州 “|| 年 平均 降水 量 Annual mean precipitation 713mm 

经 纬度 Latitude and longitude 42?24'0"N , 128?05'45"E 全 年 日 照 时 数 Annual sunshine hours 2271 一 2503h 
海拔 Altitude 738 m 全 年 无 霜 期 Annual frost free period 109— 141d 
气候 类 型 Climate type 温带 季风 气候 年 平均 气温 Annual mean temperature 3.6%C 

植被 类 型 Vegetation types its EE Bs pp WAMI Observation time 2002 年 至 今 
优势 物种 Dominant species 红 松 色 木 水 曲 柳 


1.2 模型 简介 

植被 光合 呼吸 模型 (VPRM ) 是 基于 卫 属 遥感 的 陆地 生态 系统 碳 收文 诊断 模型 , 它 由 Mahadevan 等 在 植被 
光合 模型 (VPM) 关 2 的 基础 上 发 展 而 来 所 ] 。 模 型 结构 如 图 1 所 示 ,模型 详细 介绍 见 文献 ,本 文 仅 对 该 模 
型 作 简要 介绍 。 

如 图 1 所 示 , 在 VPRM 中 ,NEE 的 计算 由 两 部 分 组 成 :由 光照 驱动 计算 的 总 生态 系统 CO, 交换 量 ( Gross 
ecosystem exchange of éarbon dioxide, GEE , Jy T Ej NEE 相对 应 ,这 里 表示 为 GEE ,其 量 值 与 GPP 相同 ,符号 相 
反 )( 式 1 等 号 右边 第 -AD) REICH SR HEP CR, ) ( 式 1 等 号 右边 第 二 、 三 项 ) 。 其 具体 表达 式 
如 下 : 


1 
x ————— 
1 + PAR/ PAR, 


AP Le ERCA F GREN HRK, PAR 为 光合 有 效 辐射 值 ( umol m? s ) ,PAR, 是 光照 为 半 饱 和 条 件 下 
光合 有 效 辆 射 值 (mol m? s!) ,FAPR, ,是 植被 吸收 的 光合 有 效 辐射 部 分 (PAV) 占 人 射 的 PAR 的 比例 ,经 研 
究 与 EVR 有 如 下 的 函数 关系 : 


NEE =- & x PAR x FAPR,y +a x T +B (1) 


FAPR,,, = EVI (2) 
模型 中 光 能 利用 率 e( pmol CO,/ pmol PPFD ) 由 最 大 光 能 利用 率 su( umol CO,/ umol PPFD) 以 及 影响 光 
能 利用 率 的 3 个 因素 ( T ue N W cale JP ) 计算 得 到 ,€ 及 T e ^ W oale .Pa 计算 公式 如 下 : 


£r £o x Teale x W sate x P ale ( 3 ) 
0 x T < Tun 或 T > T max 
Tid 7 ( T ES Trin) ( T = dum ) (4) 


当 T us < 7 < 7 


max 


(T = Tain) (T = To) Z T- P) 


min max 
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1 + LSWI 


| = T+ LSWI (5) | 卫星 遥感 资料 
1 + LSWI m 
sae ooo (6) 气象 数据 7 PAR. | EVI, LSWI 
式 中 ,7 表示 光 合作 用 的 温度 敏感 性 , WW. 表示 水 分 
对 光合 作用 的 影响 ， P. AETR AER jio 
的 影响 指数 。 方 程 (4) 中 的 Ts ,7 ,7,, 分 别 表示 光合 | 模型 参数 eo PARo a f VPRM 
作用 所 需要 的 最 低 、 最 高 和 最 适 气温 ， 单位 都 是 C , 当 气 
温 低 于 7, 时 ,7,, 取 0%C。 将 长 白山 站 2005 年 生长 季 GEE- -XI + PAR/PARY) X PAR X FAPAR pa; 
5 一 9 Hl GE AER 3C 分 组 ,分 别 计算 各 温度 段 的 CEE z; i 
和 平均 气温 ,并 通过 分 析 得 到 植物 光合 作用 最 适 气温 一 | 
7 为 28C。 另 外 ,本 研究 将 Tuns 7, 分 别 设 为 0C 和 | 
40*C , LSWI 为 地 表 水 分 指数 ,LSWI,, 为 站 点 生长 季 内 NEE = GEE} R. 
最 大 的 LSWI 值 。7 到 ， 和 已 ， 的 取 值 范围 为 0.0 一 V iiti 
1.0, NEE 


5j VPM 相 比 ,VPRM 增加 了 呼吸 项 。 将 生态 系统 
呼吸 简化 为 气温 的 线性 函数 (ax7+B)。 其 中 w (pmol 
CO, m? gc ) 和 B ( umol C 0, m? ENTE 中 为 呼吸 mer i3 diagram of the Vegetation Photosynthesis and 
Respiration Model ( VPRM) 
参数 。 这 两 个 参数 以 及 前 BU 面 的 Eo 和 PAR, 都 可 利用 模 EVI; 增 强 型 植被 指数 enhanced vegetation index;LSWI: 地 表 水 分 指 
拟 区 域 当地 的 观测 资料 进行 调整 优化 。 综 上 所 述 DCN nace water index; T; AE air temperature; PAR JE 4 DR 


Eb VPRM 模型 结构 示意 图 


VPRM 具体 表现 为 ， 辐射 photosynthetically active radiation ; so :最 大 光 能 利用 率 ( 或 最 大 
NEE=-e XT x W x EVI x 光量 子 效 率 ) the maximal light use efficiency ( the maximal light 
S 0 scale scale Pos 


quantum efficiency) ; PAR, :光照 为 半 饱 和 条 件 下 光合 有 效 辐 射 值 
xPAR+axT+ B (7) the half-saturation value of photosynthetically active radiation; œ „8 ; IÝ 
吸 参数 o, B, respiration parameters ; GEE :总 生态 统 CO; 交换 gross 


1 
1 * PAR/ PAR, 


1.3 遥感 数据 的 获取 

本 研究 所 用 的 卫星 数据 是 从 MOBDISO9A1( 空 间 分 
辩 率 为 500 m) 数据 中 提取 的 有 绪 盖 该 观测 站 点 的 数据 ， 
时 间 分 状 率 为 8 do, 使 用 蓝 色 站 (459 一 479 nm) ,红色 光 (620 一 670 nm) , 近 红 外 波段 (841 一 875 nm) 以 及 短波 
红外 波段 (1628 一 1652 nm) 这 4 个 波段 的 地 表 反 照 率 计算 得 到 增强 植被 指数 (EVI) 和 地 表 水 份 指数 (LSWI) 。 
Bars 人 SUIT; 


ecosystem CO, exchange; R, ; ^E zs & IPIK Respiration; NEE : 净 生 


态 系统 CO, 交换 net ecosystem CO, exchange 


Phir ~ Pred 
EVI =G x 2 (8) 
Pair + CC X Prea 7 Ca X Prue) +L 
Prir 7 Pswir 
LSWI = 一 一 一 (9) 
Puis + Duis 


起 中 BG C, .C, 和 工 为 常数 ,G=2.5,C,=6,C,=7.5,L=1。p 代表 相应 波段 的 地 表 反 照 率 。 下 标 nir red, blue 
和 Wwif 分 别 代 表 红 外 波段 (841 一 876 nm) 、 红 光 ( 620 一 670 nm) .蓝光 (459 一 479nm) 和 短波 红外 波段 (1628 一 
1652 nm ) 。 
1.4 通 量 观测 数据 及 小 气候 数据 的 获取 

本 研究 所 用 的 CO, 通 量 数据 是 由 开路 式 涡 度 相关 观测 得 到 的 ,该 观测 系统 位 于 样 地 内 的 气象 观测 塔 40 
m 高 度 处 ,该 系统 由 三 维 超声 风速 仪 (CSAT3 , Campbell, USA ) 和 C0,/H, O 红外 气体 分 析 仪 (Li7500 , Li-cor, 
USA) 构 成 。 观 测 从 2002 年 开始 ,数据 采样 频率 为 10 Hz, 通 量 平均 时 间 为 30 min。 观 测 塔 观测 的 详细 描述 以 
及 通 量 观测 数据 的 后 处 理 参 考 中 国 通 量 观测 研究 联盟 ( ChinaFLUX ) ( http ://www. chinaflux. org ) MRR E 
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此 外 ,气象 塔 安装 了 7 层 常 规 气 象 观测 系统 ,对 气温 、 相 对 湿度 .风速 .土壤 热 通 量 有效 辐射 .总 辐射 等 气象 要 
素 进 行 观测 。 本 文 所 需 PAR 数据 ,由 位 于 32 m 的 光 和 有 效 辐射 仪 观测 得 到 ,由 于 该 站 点 从 2002 年 开始 观 
测 ,获取 的 PAR 数据 可 能 已 经 发 生 了 衰减 。 因 此 ,实验 利用 观测 第 一 年 所 获得 数据 资料 中 的 太阳 总 辐射 与 
PAR 建立 的 关系 ,对 2006 年 PAR 数据 进行 校正 。 校 正 过 程 中 发 现 使 用 校正 前 后 的 PAR 模拟 的 NEE 数据 差 
异 并 不 显著 ,因此 ,本 文中 使 用 原始 的 PAR 观测 数据 。 
1.5 资料 预 处 理 参数 的 优化 与 验证 

参照 ChinaFLUX 以 及 李 春 等 ' 通 量 观 测 数据 处 理 方 法 和 流程 (http://www.chinaflux.org) 对 20052006 
年 长 白山 站 的 NEE 数据 进行 处 理 ,具体 步骤 包括 3 次 坐标 旋转 、WPL 校正 储存 项 计算 ,数据 质量 控制 步 又 
包括 剔除 降水 同时 段 的 数据 .剔除 夜间 摩擦 风速 低 于 临界 摩擦 风速 对 应 的 NEE 数据 .剔除 超出 帆 值 明显 异常 
的 数据 对 NEE 数据 进行 连续 5 点 3 倍 标准 差 的 剔除 。 其 中 2005 年 的 数据 用 于 模型 参数 的 优化 ;2006 年 的 
数据 用 于 VPRM 模型 优化 后 的 验证 。 

参数 优化 方案 分 为 两 步 ,第 一 步 , 利 用 2005 年 全 年 夜间 的 NEE 观测 值 通过 线性 拟 合 得 到 呼吸 参数 和 
(由 于 夜间 植物 没有 光合 作用 , 所 以 夜间 的 NEE 数据 全 部 表现 为 生态 系统 呼吸 ) 。 第 二 步 ,利用 w 和 Bp 以 及 
白天 各 时 刻 的 气温 并 计算 得 到 白天 的 生态 系统 呼吸 R., 由 白天 的 NEE 减 去 -R. 得 到 eCEE。 再 利用 VPRM 中 
GEE 的 计算 公式 ( 即 式 (7) 中 的 等 号 右边 第 一 项 ) 以 及 对 应 的 T oe WW, Poa EVI BAR 值 , 反 演 拟 合 得 到 e, 
fll PAR, 

验证 方案 分 为 3 个 部 分 。 前 两 个 部 分 :将 两 组 参数 得 到 的 NEE 模拟 结果 按 植物 生长 季 (5 一 9 H ) .植物 
非 生 长 季 (1 一 4 月 .10 一 12 H ) 两 个 的 时 间 段 分 开 , 与 观测 值 人 比较 。 统 计 分 析 方 法 包括 30 minNEE 模拟 值 
与 观测 值 的 回归 分 析 、NEE 模拟 值 与 观测 值 日 平均 变化 的 对 比 以 及 回归 分 析 。 第 三 部 分 :将 2006 年 7 月 
23—29 日 晴天 天 气 时 段 内 和 6 月 9 一 15 日 阴雨 天 气 时 段 内 NEE 的 模拟 结果 与 观测 值 作对 比分 析 。 统 计 分 析 
方法 为 30minNEE 模拟 值 与 观测 值 随时 间 变 化 的 对 比 。 


2 结果 与 分 析 


2.1 VPRM 模式 参数 优化 

利用 长 白山 站 2005 年 的 观测 资料 对 阔 叶 红 松林 的 VPRM 参数 进行 拟 合 优化 ,结果 如 表 2 所 示 。 表 2 中 
原始 参数 取 自 Mahadevan 等 在 北美 地 区 利用 Harvard. Forest 站 点 [5 温带 落叶 阔 叶 林 2000—2003 年 观测 资料 
对 VPRM 参数 的 反 演 人 %| 。 VPRWB 参 数 经 过 优化 后 ,eo 由 原先 的 0.127 变 为 0.351 ,为 原先 值 的 2.8 倍 。 同 时 ， 
PAR, 也 由 原先 的 570.0 umol m^? s^ Æ X 279.6 umol m? s^ , 减 小 到 原先 值 的 0.49 fi, B 由 原先 的 0.250 
hmolCO, m^ s’ ÆW 1.541 umolCO, ms 。 两 组 参数 存在 较 大 的 差异 主要 与 不 同 生 态 系统 植被 的 光合 生理 
特性 有 关系 ,VPRM 模式 参数 经 过 优化 后 ,更 加 能 够 反映 长 白山 地 区 的 气候 和 植被 类 型 特征 。 


表 2 VPRM 模式 参数 
Table 2 Parameters of VPRM 


参数 go/ PAR,/ a/ " 
eee -2 .-1q^-1 B/ (umol m s!) 

Parameters (HmolCO,/ kmolPPFD ) ( umol m? s^!) (hmolCO m^s C7) 

原始 Origunal 0.127 570.0 0.271 0.250 

优 货 后 Optimized 0.351 279.6 0.246 1.541 


2.2 VPRM 模拟 结果 与 分 析 
由 于 长 白山 地 区 森林 NEE 在 生长 季 与 非 生长 季 表 现 出 明显 不 同 的 变化 特征 ,将 2006 年 NEE 的 模拟 结 
果 分 为 植物 生长 季 (5 一 9 月 ) 与 植物 非 生长 季 (1 一 4 月 .10 一 12 月 ) 两 个 时 段 与 观测 值 进行 对 比 。 
2.2.1 植物 生长 季 NEE 模拟 值 与 观测 值 的 对 比 
将 两 组 参数 输入 VPRM 模式 ,对 2006 年 长 白山 地 区 阔 叶 红 松林 30minNEE 进行 模拟 。 对 比 5 至 9 月 ( 植 
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物 生 长 季 )30minNEE 的 模拟 值 与 观测 值 可 以 看 出 (图 2) ,在 植物 生长 季 ,参数 优化 后 模拟 的 NEE 比 采 用 原始 
参数 模拟 出 的 NEE 与 观测 值 更 为 吻合 。 原 始 参数 低估 了 白天 和 夜晚 的 NEE 值 ,回归 方程 斜率 仅 为 0.488 ,在 
NEE 值 大 于 10 pmol m? s^ 和 小 于 -20 pmol m^ s^! 时 ,低估 尤为 明显 。 参 数 优化 后 ,模型 对 NEE 的 低估 得 到 
改善 ,回归 方程 斜率 达到 0.911。 采 用 原始 参数 模拟 的 NEE 值 与 观测 值 的 平均 误差 为 3.19 umol m^ s^ , 均 方 
根 误差 为 6.76 umol m? s! ,参数 优化 后 NEE 模拟 值 与 观测 值 的 平均 误差 为 -1.81 pmol m^? s^! , 均 方 根 误差 
为 8.61 pmol m^ s (2€ 3) 。 人 参数 优化 后 模拟 值 与 观测 值 的 均 方 根 误差 稍稍 变 大 ,但 平均 误差 减 小 ,模拟 结 
总 体 优 于 参数 优化 之 前 。 


30 p 30 p 
a b 
mc 20 F E 20r 
RS D RS o 
SÉ g 10 上 ru 10 - 
X3 0 X o 
$S -10 y SE -10 } 
= 0.488x + 0.62 A " 
$Z -20 J= 0488x0620, ut 20 h= 0 leet 253 
Re Ra ~ 
> > E 
Mu 99 m -30 
Ka d 
xb -40 HB -40 
—-50 1 1 1 L 1 1 1 -50 L 1 1 1 1 1 1 1 
-50 -40 -30 -20 -10 0 10 20 30 -50 -40 -30 -20 -10 0 10 20 30 
EDR BER CO» AE BORN ER 


NEE obs/(umol m? s!) 


2 VPRM 原始 参数 (a) 及 优化 后 参数 (b ) 模拟 的 长 白山 站 2006 年 植物 生长 季 (5S 一 9 H ) 30minNEE ( NEE. VPRM) 与 观测 值 ( NEE_ 
Obs) 的 回归 分 析 
Fig.2 Regression analyses between VPRM-simulated NEE ( NEE VPRM) and observed NEE ( NEE Obs) during growing season ( May 


to September) in 2006 at Changbaishan site based on original parameters (a) and optimized parameters (b) 


对 生长 季 NEE 模拟 值 与 观测 值 的 平均 日 变化 进行 "AN 
对 比 可 以 看 出 (图 3) , VPRM 较 好 地 模拟 出 也 INEE 单 xl iol 

峰 型 的 变化 趋势 ,但 对 NEE 峰值 存在 一 定 前 低估 。 原 | m "m 
始 参 数 模拟 的 NEE 峰值 为 -6.32 iumol. m? s, RÆ kp Mea = Ka 
测 值 峰值 ( -14.66 umol m? s~!) 的 3%\ 优 化 后 参数 模 E i T 2 

拟 得 到 的 NEE 峰值 为 -13.27 umo m? s^, WAWE EF 5} e ML mut 

d 9190, IG FEL AE RCM 相 比 优化 之 前 , 模 Sus AMA 

拟 值 与 观测 值 的 平均 误差 减 小 到 -0.06 pmol m? 57,38. — | ro 

方 根 误 差 下 降 到 2.17Wumol m? s! PERKER 3), EP SEP SEP SEP SEP SED SEE SE SE SEP SEE MEUS 
模拟 得 到 明显 改善 。 另 外 ,模型 对 夜间 森林 的 呼吸 存在 三 
— 5E BA If.) 主要 原因 与 模型 使 用 线性 的 方程 估计 R, TER 

有 很 大 的 关系 。 图 3 长 白山 站 2006 年 植物 生长 季 (5 一 9) 月 模拟 的 NEE 与 观测 


Bb PKE NEE 平均 日 变化 模拟 值 与 观测 值 的 平均 日 变化 

值 的 回归 分 析 图 ( 图 4) A 可 以 看 出 ， 采 用 优 1 化 后 E 数 模 Fig.3 A comparison between the observed and simulated mean 
diurnal variation of NEE during growing season ( May to 

e d ibd Bd n A 附近 ° 相 较 之 下 i 采用 原始 参 September) in 2006 at Changbaishan site 
数 得 到 的 模拟 值 与 观测 值 的 回归 方程 斜率 仅 为 0.497， NEE Obs. NEE I observed NEE; NEE. Ori. 原始 参数 模拟 的 
模拟 效果 相对 较 差 。 该 结果 进一步 表明 ,优化 后 的 NEE simulated NEE with original parameters;NEE_Opt: 优 化 后 参数 
VPRM 能 更 好 地 模拟 植物 生长 季 NEE 的 日 变化 。 模拟 的 NEE simulated NEE with optimized parameters 
2.2.2 ”植物 非 生长 季 NEE 模拟 值 与 观测 值 的 对 比 

图 $ 为 植物 非 生 长 季 (1 一 4 月 10 一 12 月 )30 minNEE 的 模拟 值 与 观测 值 的 对 比 。 由 图 可 以 看 出 使 用 两 
组 参数 VPRM 模拟 非 生长 季 30 minNEE 的 效果 都 不 理想 。 原 始 参数 和 优化 后 参数 各 自 对 应 模拟 结果 与 观测 
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图 4 长 白山 站 2006 年 植物 生长 季 (5 一 9 月 )NEE 的 模拟 值 (NEE_VPRM ) 与 观测 值 ( NEE. Obs) 平均 日 变化 的 回归 分 析 
Fig.4 Regression analyses between the simulated mean diurnal variation of NEE ( NEE-VPRM) and observed mean diurnal variation of 


NEE (NEE-Obs) during growing season ( May to September) in 2006 at Changbaishan site 
NEE Ori, 原始 参数 模拟 的 NEE simulated NEE with original. parameters; NEE |. Opt, 优化 后 参数 模拟 的 NEE simulated. NEE with 


optimized parameters 


值 拟 合 方程 的 斜率 分 别 为 0.058 和 0.068 ,模拟 对 NEE 高 箱 都 存在 明显 的 低估 。 采 用 原始 参数 和 优化 后 参数 
模拟 的 NEE 值 与 观测 值 的 平均 误差 分 别 为 -0%86\WumoIn- s 和 0.39 umol m? s! , RIRIA aA 0.005 和 
0.01( 表 4) ,模拟 效果 都 较 差 ,这 主要 是 因为 在 植物 非 生 长 季 气 温 比 较 低 ,此 时 的 生态 系统 呼吸 与 气温 之 间 不 
符合 线性 关系 1 ,而 VPRM 把 呼吸 简化 为 气温 的 线性 函数 ,这 将 导致 模型 在 植物 非 生长 季 对 NEE 的 模拟 
存在 较 大 的 误差 。 


表 3. 2006 年 植物 生长 季 ( 5 一 9 月 ) NEE 模拟 值 与 观测 值 之 间 的 统计 分 析 
Table 3 Statistical analysis between simulated NEE and observed NEE during plant growing season ( May to September) in 2006 


均 方 根 误差 


数据 类 型 采用 两 组 不 同 的 参数 平均 误差 Fin 相关 系数 的 平 有 效 数 据 
Data types Using two different parameters Average error Bud rror FIR R? Valid data 
30minNEE 采用 原始 参数 3.19 6.76 0.53 4115 
采用 优化 后 参数 -1.81 8.61 0.52 4115 
NEE 日 平均 变化 采用 原始 参数 2.65 4.31 0.93 48 
Mean diurnal variation-of NEE 采用 优化 后 参数 —0.06 2.17 0.93 48 


RA 2006 年 植物 非 生长 季 NEE 模拟 值 与 观测 值 之 间 的 统计 分 析 

Table 4 In 2006, statistical analysis between simulated NEE and observed NEE during non-growing season ( January to April and October to 
December) in year 2006 

均 方 根 误差 


数据 类 型 采用 两 组 不 同 的 参数 平均 误差 Rome 相关 系数 的 平 有 效 数 据 
Data types Using two different parameters Average error square emol 方 值 R? Valid data 
30minNEE 采用 原始 参数 -0.86 3.34 0.005 5372 
采用 优化 后 参数 0.39 3.09 0.01 5372 
NEE 日 平均 变化 采用 原始 参数 -1.24 0.98 0.001 48 
Mean diurnal variation of NEE 采用 优化 后 参数 0.01 0.80 0.25 48 
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5 VPRM 模型 基于 原始 参数 (a) 与 优化 参数 (b) 模拟 的 长 白山 站 2006 年 非 植物 生长 季 ( 1 一 4 H 10-12 H ) 80minNEE(NEE. VPRM) 
与 观测 值 ( NEE_Obs) 的 回归 分 析 
Fig.5 Regression analyses between VPRM \simulated NEE (NEE_VPRM) and observed NEE ( NEE-Obs) during non-growing season 


(January to April and October to December) in year 2006 at Changbaishan site based on default parameters (a) and optimized parameters 


(b) 


比较 植物 非 生长 季 NEE 平均 日 变化 的 模拟 值 与 观测 值 发 现 (图 的: 植物 非 生 长 季 NEE 的 观测 峰值 ( 约 
为 -1.70 umol ms ) 比 生长 季 要 小 的 多 , 仅 为 生长 季 峰 值 的 12965 Je ARCU] NEE 峰值 为 -1.60 umol 
ms ,与 观测 值 较为 接近 ,但 模拟 在 上 午 和 下 午 都 存在 严重 的 高 和 估 。 参 数 优化 后 模拟 得 到 的 NEE 峰值 仅 为 
-0.300 pmol ms”, 相 比 观测 值 低估 了 82%。 同 时 植物 非 生长 季 夜 间 呼 吸 值 (平均 为 0.993 kmol m^ s) th, 
低 于 生长 季 ( 平 均 为 5.77 pmol m”? s) 。 原 始 参 数 模拟 夜间 呼吸 值 为 负 值 ,其 模拟 结果 不 可 信 。 优 化 参数 后 
模拟 夜间 NEE 的 平均 值 为 0.367 umol m^ s^! , 相 比 观测 值 低 估 了 63% 。 在 整体 上 ,参数 优化 后 模拟 非 生长 季 
NEE 平均 日 变化 的 平均 误差 为 0.01 umol m? s^ IN HR E25 N 0.80 pmol m? s^! ,与 观测 数据 的 R7 H 0.25, 
而 原始 参数 的 模拟 效果 在 这 3 项 指标 上 均 不 如 优化 后 ( 表 4) 。 

对 植物 非 生长 季 NEE 日 平均 变化 模拟 值 与 观测 值 进行 回归 分 析 可 以 看 出 (图 7) ,两 组 参数 对 植物 非 生 
长 季 NEE 的 模拟 效果 都 较 差 。 使 用 原始 参数 和 优化 后 参数 模拟 的 NEE 值 与 观测 值 回归 方程 斜率 仅 为 0.012 
和 0.175。 以 上 分 析 表 明 , 虽 然 参 数 优化 后 的 模拟 方案 优 于 优化 之 前 ,但 是 VPRM 模型 对 植物 非 生长 季 NEE 
的 模拟 效果 较 差 。 

从 以 上 不 同 生 长 肥 NEB 模 拟 值 与 观测 值 的 对 比 可 表明 , 相 比 原始 参数 ,由 观测 数据 拟 合 的 参数 更 加 适用 
于 VPRM 对 长 白山 地 区 NEE 的 模拟 ,尤其 是 生长 季节 。 因 此 ,文章 后 面 只 讨论 使 用 优化 参数 模拟 得 到 的 模拟 
值 与 观测 值 的 对 比 。 

2.2.3 ”阴雨 天 握 对 模拟 的 影响 

由 ,NEB 观 测 数据 可 知 ,阴雨 天 气 时 由 于 云层 的 阻挡 ,白天 入 射 到 地 面 的 太阳 辐射 会 降低 ,这 会 导致 NEE 
减 小 ,而 在 夜间 过 层 存在 向 下 逆 辐 射 会 导致 地 面 增 温 ,这 可 能 会 导致 呼吸 作用 的 增加 。 因 此 ,有 必要 分 析 阴 两 
天 气 对 IENFEE 模拟 的 影响 。 

将 2006 年 7 月 23—29 H (DOY204-210) 晴天 天 气 和 6 月 9—15 日 (DOY160- 166) 阴雨 天 气 VPRM 的 模 
拟 结果 与 观测 值 作对 比 ( 图 8)( 根 据 日 最 大 光 和 有 效 辐 射 值 划 分 晴天 和 阴雨 天 :晴天 >1200 umol m^ s~, BH 
雨天 <1200 umol m° s); HUA H, ERRAT (E 8a), VPRM 较 好 地 模拟 出 了 NEE 的 变化 动态 ， 
NEE 观测 值 在 白天 的 峰值 为 -18.32 一 -29.48umol ms ,模拟 值 的 峰值 为 -27.52 一 -33.95 umol m^ s^! ,模拟 
对 NEE 的 峰值 略 有 高 估 , 这 是 因为 模型 并 没有 考虑 ,晴天 正午 太阳 辐射 强度 较 高 的 条 件 下 ,导致 的 植被 “光合 
午睡 ”的 现象 。 而 在 阴雨 天 条 件 下 (图 8b) , NEE 观测 值 波动 较 大 ,NEE 白天 峰值 为 -18.45 一 -24.67hmol m^ 
s ,和 白天 模拟 值 峰 值 为 -4.77 一 -10.14umol m? s, VPRM 不 仅 没 有 模拟 出 NEE 的 变化 趋势 ,而 且 对 白天 峰 
值 存在 明显 的 低估 ,模拟 效果 较 差 。 这 表明 VPRM 模型 对 阴雨 天 NEE 的 模拟 还 需要 进一步 改进 。 
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图 7 长 白山 站 2006 FEHI 2k 1€ 28 (1—4 H..10—12 月 )NEE 的 
模拟 值 ( NEE_VPRM) E3)HBMÉ CNEE Obs) 325 A S£ 4£ B [903 23 f 
Fig.7 Regression analyses between the simulated mean diurnal 


variation of NEE ( NEE-VPRM ) and observed NEE ( NEE-Obs ) 


6 长 白山 站 2006 年 植物 非 生长 季 模 拟 的 NEE 与 观测 值 的 平 

均 日 变化 

Fig.6 A comparison between the observed and simulated mean 

diurnal variation of NEE during non-growing season of 2006 at 
. : during non-growing season ( January to April and October to 

Changbaishan site 

NEE, Obs, NEE 观测 值 observed NEE; NEE, Ori, 原始 参数 模拟 的 

NEE simulated NEE with original parameters; NEE, Opt , 优化 后 参数 


December) of 2006 at Changbaishan site 
NEE Om, 原 始 参 数 模 拟 的 NEE simulated NEE with original 


" arameters; NEE. Opt, 优化 后 参数 模拟 的 NEE simulated. NEE with 
模拟 的 NEE simulated NEE with optimized parameters para JAg zOPt, 优化 后 参数 模拟 的 e m 


optimized parameters 


图 9 为 晴天 和 阴雨 天 VPRM 模拟 的 GEE 与 尺 . 的 变化 值 。 晴 天 GEE 的 峰值 逐日 上 升 ,R. 的 变化 保持 稳 
定 , 日 变化 幅度 较 GEE 小 。GEE 的 峰值 为 -33.5 针 ~-41.83 umol m° s^! ,R 为 4.05 一 7.94umol m° s~, HERF, 
NEE 主要 受 GEE 控制 ,这 体现 了 选取 合适 的 俩 和 PAR, 值 对 晴天 条 件 下 模拟 NEE 的 重要 性 。 同 时 可 看 出 R, 
的 峰值 相 比 GEE 的 峰值 延迟 约 1 一 2h ,这 十 因 擅 WR. 主要 受气 温 控 制 , 近 地 面 气温 在 午后 2h 左右 达到 最 高 ,而 
植物 光合 作用 所 需 的 太阳 辐射 最 强 的 时 间 症 在 正午 。VPRM 模型 较 好 地 模拟 了 上 述 过 程 。 
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8 长 白山 站 2006 ££ DOY204—210( 晴天 ) (a) E DOY160—166( 阴雨 天 )(b)NEE 模拟 值 与 观测 值 的 对 比 
Fig.8 A comparison between simulated NEE and observed NEE at Changbaishan site in DOY204-210 (sunny days) (a) and DOYI160-166 
(cloudy or rainy days) (b) in 2006 
NEE-Obs , NEE 观测 值 observed NEE ; NEE. VPRM ,参数 优化 后 模拟 的 NEE simulated NEE with optimized parameters 


与 晴天 相 比 ,阴雨 天 GEE 值 较 小 旦 变化 起 伏 较 大 ,R. 变 化 无 规律 (图 9b)。GEE 的 峰值 为 -9.96 一 -16.89 
umol ms ,RR. 为 3.77 一 6.98 umol m^ s! , GEE 的 最 大 峰值 与 R. 最 大 值 都 出 现在 第 163 天 ,当日 PAR 的 最 
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大 值 为 1632 umol m? s ,而 GEE 最 小 峰值 出 现在 第 162 天 ,当日 PAR 峰值 为 622 umol m° s, E FHR 
气 的 影响 ,云层 对 太阳 辐射 的 遮挡 作用 ,PAR 在 模拟 时 间 内 变化 起 伏 大 ,导致 GEE 变化 起 伏 大 。R. 变 化 不 规 
律 是 因为 生态 系统 呼吸 主要 由 气温 (或 土壤 温度 ) 控 制 ,而 在 阴雨 天 时 ,云层 白天 对 太阳 辐射 有 削弱 和 夜晚 向 
下 逆 辐 射 对 地 面 存在 增 温 ,地 面 气温 (或 土壤 温度 ) 变化 日 较 差 减 小 ,从 而 导致 生态 系统 呼吸 相应 的 变化 。 
此 ,VPRM 对 阴雨 天 NEE 模拟 存在 较 大 的 不 确定 性 。 

综 上 所 述 , 通 过 对 不 同 天 气 个 例 的 分 析 , 发 现 VPRM 模型 可 以 较 好 的 模拟 晴天 天 气 下 长 白山 地 区 NEBE 的 
变化 ,但 是 模型 对 阴雨 天 NEE 的 模拟 仍 显 不 足 ,模拟 的 不 确定 性 较 大 。 在 2006 年 ,长 白山 地 区 植物 生 芍 季 阴 
雨天 天 数 长 达 63 d, 占 生长 季 总 天 数 的 比率 41.2%, VPRM 模型 对 阴雨 天 NEE 模拟 的 不 足 会 影响 生长 季 
NEE 的 模拟 效果 。 而 植物 生长 季 正 是 该 地 区 森林 NEE 发 生 的 主要 时 期 ,这 将 进一步 影响 VPRM 对 全 年 NEE 
的 模拟 。 
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9 长 白山 站 2006 年 DOY204- 210( 上 晴天) (a) 与 DOY160-166( 阴雨 天 )(b) 的 GEE 5 RK. 的 变化 
Fig.9 A comparison between simulated and observed GEE and R, at Changbaishan site during DOY204-210 (sunny days) (a) and during 
DOYI160-166 (cloudy or rainy days) (b) in 2006 
GEE ,参数 优化 后 模拟 的 GEE simulated GEE with optimized parameters ; R, ,参数 优化 后 模拟 的 R, simulated respiration with optimized parameters 


3 讨论 


最 大 光 能 利用 率 so 是 VPRM 估计 GEE 大 小 的 关键 参数 ,而 不 同 植被 类 型 或 气候 条 件 的 区 域 所 对 应 的 es 
并 不 同 。 本 研究 利用 长 白山 站 2005 年 通 量 观测 资料 反 演 拟 合 得 到 的 温带 阔 叶 红 松林 的 so 值 为 0.351。 比 较 
类 似 研 究 的 s, 值 ;例如 Mahadevan 等 反 演 北美 HARVARD 站 点 ' 引 观测 资料 所 得 到 的 e, fti 7g 0.1277! , Yuan 
等 得 到 的 VPRM 落叶 阅 叶 林 的 e EA 0.156! , Hilton 等 得 到 的 落叶 阔 叶 林 se, 平 均值 为 0.24!:29 ,本 研究 得 
到 的 su 值 更 大 这 可 能 与 长 白山 阔 叶 红 松林 位 于 季风 区 ,生长 季 内 , 光 ,温水 充足 这 会 有 利于 该 区 域 植被 的 
光合 作用 ,着 提高 其 光 能 利用 率 。 

参数 优化 之 后 ,VPRM 可 以 较 好 地 模拟 植物 生长 季 NEE 的 变化 ,但 对 植物 非 生 长 季 NEE 的 模拟 效果 较 
六 $ 这 与 刘 诚 等 人 的 研究 结果 相同 :2 。 造 成 这 样 差异 的 原因 是 模型 把 生态 系统 呼吸 仅 作为 气温 的 线性 函数 。 
目前 ,生态 系统 呼吸 随 着 温度 呈 指 数 增加 趋势 的 认识 得 到 了 广泛 认同 。 同 时 也 有 研究 表明 ,温带 地 区 森林 生 
态 系统 呼吸 与 土壤 表层 温度 的 相关 性 最 强 , 可 以 解释 生态 系统 呼吸 变异 的 85% 左 右 , 而 近 地 面 温度 只 能 解释 
生态 系统 呼吸 变异 约 70%5 。 土 壤 温 度 与 生态 系统 呼吸 有 着 更 好 的 关系 。 因 此 在 未 来 ,可 以 针对 温带 森林 
生态 系统 对 模型 进行 改进 ,用 指数 关系 代替 原本 的 线性 关系 ,用 土壤 温度 代替 气温 来 模拟 生态 系统 呼吸 从 而 
改善 VPRM 模型 对 NEE 的 模拟 。 另 外 ,在 阴雨 天 条 件 下 ,模型 的 模拟 也 存在 一 定 的 不 确定 性 。 阴 雨天 气 条 件 
下 ,土壤 售 水 量 的 变化 对 生态 系统 呼吸 的 影响 以 及 生态 系统 散射 辐射 的 增强 导致 植物 的 光合 固 碳 过 程 的 变化 
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VPRM 并 没有 考虑 ,这 也 是 在 未 来 模型 结构 需要 改进 的 地 方 。 

除了 模型 结构 的 不 足 ,在 实验 数据 的 获取 过 程 中 ,也 存在 一 定 的 不 确定 性 。 由 于 阴雨 天 时 受 云 层 的 遮挡 
等 因素 影响 ,获取 的 遥感 数据 会 存在 一 定 的 噪音 :5 。 这 些 都 会 导致 模型 对 阴雨 天 NEE 的 模拟 存在 不 确 
定性 。 

在 其 他 VPRM 模型 参数 优化 及 验证 的 研究 中 ,如 Mahadevan 45 X Hilton 等 人 分 别 使 用 了 北美 地 区 22 .65 
个 站 点 2 一 7a 的 观测 资料 ,参数 优化 后 不 同 植被 类 型 VPRM 的 模拟 效果 各 有 差异 ,个 别 站 点 模拟 效果 不 理 
48/77". Mahadevan 等 人 使 用 优化 后 的 VPRM 模型 模拟 HARVARD 站 点 温带 阔 叶 林 森林 生长 季 NEE/ 的 误差 
为 -0.26 pmol m^ s^! ,相关 系数 为 0.91 ,结果 优 于 本 研究 。 不 过 他 们 的 研究 没有 涉及 到 非 生长 季 NEE BU, 
阴雨 天 对 模拟 的 影响 及 R. 与 温度 的 线性 关系 问题 。 本 研究 只 使 用 了 长 白山 一 个 站 点 两 年 的 观测 资料 ,人 研究 
若 能 获得 更 多 站 点 更 长 时 间 的 观测 资料 ,参数 的 优化 及 结果 能 够 体现 出 气候 及 植被 生长 的 变 北 , 将 有 助 于 改 
善 VPRM 对 个 别 站 点 NEE 模拟 的 不 足 。 我 国 幅 员 辽 阔 , 气 候 和 植被 类 型 多 样 ,大 能 将 VPRMI 应 用 到 更 多 地 
域 并 进行 进一步 的 优化 ,得 到 一 组 适合 我 国 不 同 植被 类 型 的 模型 参数 ,这 对 理解 我 国 区 域 陆 屯 矶 循环 过 程 具 
有 重要 意义 。 


4 结论 


本 研究 利用 ChinaFLUX 长 白山 站 阔 叶 红 松林 2005 年 通 量 观测 资料 对 VPRM 模型 关键 参数 (se PAR, 、a 
和 6) 进行 了 优化 ,通过 对 阔 叶 红 松林 2006 年 植物 生长 季 (5 一 9 月 ) 和 植物 非 生长 季 (1 一 4 月 与 10 一 12 H) 
NEE 的 模拟 ,对 优化 得 到 的 参数 进行 了 验证 ,并 分 析 了 不 同 天 氢 条 件 对 /NEE 模拟 的 影响 ,得 到 以 下 结论 : 

(1) 对 于 温带 针 阔 叶 混 交 林 , VPRM 模型 中 的 e, PAR wa 4 8 经 优化 后 的 值 为 0.351 pmolC0,/ 
umolPPFD ,279.6 umol m ? s™' .0.246 umolCO, m ? s ' ‘CHA 和 1.541 pmol m? s^! , 

(2) 优 化 后 的 VPRM 能 够 较 好 地 模拟 长 白山 地 区 植物 生长 季 NEE 的 变化 。 对 30minNEE 模拟 的 平均 误 
差 为 -1.81 pmol m^ s^ ,相关 系数 为 0.72。VPRM 模 型 对 植物 非 生长 季 森 林 NEE 的 模拟 效果 较 差 。 模 拟 
30minNEE 的 平均 误差 为 0.39 pmol ms 与 相 关系 数 仅 为 0.10 ,并 且 模 拟 低估 NEE 平均 日 变化 白天 吸收 峰值 
约 82% ,日 变化 模拟 值 与 观测 值 的 相关 系数 仅 为 0.50。 

(3) 通 过 对 不 同 天 气 个 例 的 分 析 , 发现 VPRM 可 以 较 好 地 模拟 晴天 条 件 下 长 白山 地 区 NEE 变化 ,但 对 阴 
雨天 NEE 的 模拟 仍 显 不 足 。 
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